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Introduction 


In  a previous  report  [l],  a complete  theory  was  presented 
of  the  propagation  of  weak  spherically  symmetric  shocks  in  air 
of  uniform  density,  pressure  etc:  Later,  the  theory  was  extended 


in  which  the  undisturbed  fluid  into  'which 


shock 


moved  was  not  uniform,  but  the  property  of  spherical  symmetry 
was  preserved.  In  the  present  report,  a combination  of  the  meth- 
ods of  [l]  and  [2]  is  applied  to  the  problem  of  a shock  (produced 
by  an  explosion  say)  propagating  in  the.  atmosphere,  taking  into 
account  the  variation  of  the  density,  pressure'  etc.  of  the  air 
with  height  above  the  ground.  Since  the  air  is  stratified  and, 
in  the  most  important  problem,  the  shock  is  initially  spherical 
(approximately)  only  the  symmetry  about  a vertical  axis  remains ; 
the  physical  quantities  are  now  functions  of  two  space  variables 
and  the  time.  The  introduction  of  an  additional  independent  var- 


iable into  the  problem  is  far  from  a trivial  extension,  and  if 
the  theory  can  be  worked  out  successfully  in  this  particular 
case,  4-t  will  st.rvc  as  a prototype  for  many  other  investigations 
(in  the  theory  of  three  dimensional  steady  supersonic  flow,  for 
example) . The  present  account  constitutes  a "progress  report" 
since  certain  points  require  further  consideration  and  justifica- 
tion before  a final  report  can  be  written. 

Only  weak  disturbances  art  treated  here,  and  as  before  the 
linearized  (acoustic)  theory  forms  the  starting  point  of  the  method. 
The  linear  results  are  inadequate  in  certain  ways,  as  explained  in 
Cl],  but  their  shortcomings  can  be  remedied  to  give  a valid  descrip- 
tion of  the  flow.  The  work,  therefore,  falls  into  two  sections;  (I) 
the  d< rivation  and  consideration  of  the  appropriate  results  of  linear 
tneory,  and  (II)  the  improvement  of  those  results  and  the  determina- 


tion of  the  shoe:..  The  first  section,  strangely  enough,  is  the  nore 
troublesome  and  the  more  incomplete , since  explicit  solutions  of  the 
linear  equations  are  not  known  in  general,  and.  wo  have  to  bo  satis- 
fied with  a certain  expansion  solution.  This  expansion 


» 
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however,  is  of  considerable  interest  for  its  own  sake  and  is  of 
general  application  to  linear  hyperbolic  equations  of  second  order. 
The  first  section  extends  the  methods  of  treating  the  linear  theory, 
which  were  developed  in  [2],  to  include  an  additional  space  variable 
(the  treatment  of  linear  theory  is  trivial  in  the  problem  of  [l], 
since  the  explicit  solution  of  spherical  sound  waves  is  well-known), 
while  the  second  section  extends  the  applicability  of  the  basic 
ideas  described  in  [lj, 

I«  The  linearized  theory 

'The  physical  quantities  are  functions  of  t\v  distance,  r,  from 
the  (vertical)  axis  of  symmetry,  the  height,  z,  above  the  ground, 
and  the  time  t.  If  u and  w are  the  velocity  components  in  the  di- 
rections of  the  r and  z axes,  respectively,  the  density  is  p and 
the  pressure  is  p,  then  the  equations  of  momentum  (including  the  ac- 
celeration of  gravity  g)  are 

Ut  + UUr  + WUZ  = -i  pr 

wt  + uwr  + = -p  pz”6' 

and  the  equation  of  continuity  is 

(3)  Pt  + upr  + v/pz  + p(ur  + p + wz)  = 0. 

It  may  be  assumed,  since  entropy  changes  are  negligible  for  weak 
disturbances,  that  p is  a given  function  of  p*  Then,  introducing 
the  velocity  potential  >#  equations  (l)  and  (2)  may  be  integrated 
to  give  Bernoulli's  equation 

(ip)  'j>t  + |'(4)p  + 4>z)  + Sz  = constant. 

In  equilibrium  u = v = 0,  and  p = F(z),  p = R(z),  a = / dp/dp  = 

A( z ) , say, 

and  from  (2)  this  equilibrium  configuration  is  determined  by 
(5)  P’ (z)  + gR{ z)  = 0. 

To  obtain  the  linearized  equation  of  motion,  th,.  -1  evi.il  a 'ms  >f  j> , p , 
p,  and  a from  their  equilibrium  values  (the  undisturbed  value 
of  »J>  is  z.ro)  are  assumed  to  be  small,  and  In  the  above  equations 
only  the  first  order  terms  in  these  small  quantities  are  retained. 


(1) 

(2) 
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Equation  (U)  becomes 


- P ( z ) = -R(z)  <t>  t. 


and  (3)  be cones 


(7)  + A +^+i-  ^-(R<j>  ) = 0. 

' f ' R frr  r R Tz' 

p ? 

In  addition,  since  a^  = dp/dp,  p.  (p-P) +/a  ( z ) j hence,  from 

(6),  (7)  is 

(8)  K.+-.  = a2^z>  3 + F^r  + ♦**  + R^V"  >• 


. Rf(z) 
’zz  + R TzT 


The  solution  of  this  equation  for  ^ is  not  known,  'out  an  ex- 
pansion valid  near  the  front  „.f  the  disturbance  will  prove  to 
be  of  value.  To  motivate  this  expansion  we  consider  two  special 
cases:  (i)  the  case  of  a uniform  atmosphere  with  i\  equal  to  a con- 
stant A and  4>  a function  of  r alone  (this  is  the  problem  of  cy*. 
o 

lindrical  sound  waves),  and  (ii)  the  case  of  plane  waves  propagat- 
ing vertically  so  that  ^ is  a function  only  of  z.  In  each  of 
these  problems  only  one  space  variable  appears,  and  the  equation 
reduces  to  the  type  considered  in  [R]j  ( i)  and  (ii)  are  now  con- 
sidered in  turn. 

( i ) Cylindrical  sound  waves 

In  this  case  (8)  reduces  to 


4,  + in  = 

frr  r ,-c  * 


and  it  is  w^  11 -known  that  the  solution  which  represents  an  out- 
going wave  starting  at  t = 0 is 

(10)  <t>  = K (pr/AQ)  h(t)  = 

Whore  h(t)  is  an  arbitrary  function,  K is  the  Bessel  function,  and 
p don ’t-.s  the  operator  in  the  Heaviside  calculus.  (T.n  .res sure 
has  also  been  denoted  by  p,  but  no  confusion  will  ensue,  since  the 
pressure  does  not  appear  again  in  this  section) . 

It  is  f ound  in  those  pr  ;bl<..:aa  (see  [l]  and  [<2] ) that  in  order 
to  determine  the  head  shock  the  linearized  solution  is  required  in 
the  region  where  A£/r  is  small ; here,  £,  is  the  characteristic  vari- 


able  t-r/A  vihich  measures  the  time  after  the  arrival  of  the  first 

0 

disturbance  (which  travels  on  4 = 0)  . Although  small  4 corresponds 
to  the  front  of  the  disturbance  in  the  linearized  theory,  the  solu- 
tion for  small  4 is  not  sufficient.  For  in  the  improved  (non-lin- 
ear)  theory,  characteristics  are  continually  meeting  the  shock  so 
that  eventually  at  large  distances  the  appropriate  value  of  the 
characteristic  variable  near  the  shock  is  no  longer  small.  Hence, 
roughly  speaking,  we  require  the  solution  for  small  4 and,  in  addi- 
tion, the  solution  for  large  r$  more  precisely,  the  solutions  for 
small  d,£,/r  is  needed. 

The  expansion  can  be  obtained  from  the  integral  in  (10),  but 
is  found  more  quickly  from  the  Heaviside  representation  of  the  so- 
lution. Small  values  of  4 correspond  to  large  values  of  p in  the 
Heaviside  representation j in  fact,  l/p  plays  a role  corresponding 
to  4 so  that  the  expansion  of  cf>  for  small  values  of  A..4/r  is  foun-~- 
f r :>m  the  Heaviside  representation  by  assuming  A /pr  is  small.  Ex- 
panding the  Bessel  function  for  large  values  of  its  argument, 
pr/A  , and  interpreting  term  by  terra,  we  have 


(11)  k 


/2a.  e-P^o  j 


J 1 1 Ao 

^ 1 ” H pr  12 


h(t) 


V 2pr  - - I - 

A - a 1 

= at-r/Ao)  - J(di)2  f1(t-r/Ao)  + jf^A)2  — , 


who  re 


f(t)  = 


an  - i*  ( t)  dxn->tus  the  nth  rejW.tx  d integral  from  0 to  t.  To  a 
first  .r.ler  approximation  .nly  the  first  term  in  (11)  noc  ' be  re- 
tained} using  it,  the  linearize  tVu  ry  may  be  imp-r  >\n  1 and  the 
head  sh..ck  -determined  by  the  pr  ce.'.urc  described  in  ( 1 ] . The  de- 
tails arc  not  given  here,  since  this  special  case  is  nly  use  1 to 
m tivato;  the  later  treatment  of  the  linear  equation  (6);  in  any 
case,  the  results  are  very  similar  t these  of  [ 3 J . 

(ii)  Vertical  pr  -go-gat Ion  of  p.i.  vie  waves 


It  is  assumed  n >w  that  the  a.ir  is  poly  tropic  i.<_.  the  . res sure 
is  aqu  ;1  t ■ kp'g  where  k and  y are  constants.  It  follows  th.  n from 


(5)  that  R(z)  cc 


(zn-zj^r, 


-h- 

md  therefore  that 


(12) 


A = A1(zo-z)2, 


whore  z„  and  An  arc  constant";,  ■‘is  z 

o 1 


zQ,  R and  A tend  to  zero,  but 


of  course  the  simple  polytrooic  ceases  to  apply  to  the  atmosphere  bc- 

2 

fere  z — zQ  is  reached  (for  example,  the  tropopause,  at  which  T cc  4 
is  discontinuous,  intervenes)}  however,  the  behaviour  according  to  (12) 
of  A near  z = zq  will  prove  interesting  mathematically.  With  these 
values  of  R and  A,  and  A a function  of  z alone,  (8)  reduces  to 


(13) 


*tt  = Al(zo“z) 


It  is  convenient  to  introduce 
(H4-) 


i Az  1 

y-i  z0-z  r 

2AI1(zo^)^ 


in  place  of  z;  then 


(15) 


1 _ * 2yAw 

\t  - y>i“  -yt  * 


w'lorc  v = ^(  3-y)/(  Y“l) » and  the  characteristics  are  t + ^ = constant. 
Replacing  A^t  by  the  operational  form  p^A,  the  general  solution  of  (15)  i 


(16) 


4 (THr  ^p,hi(t)  + 


whore  I.. , 1 and  K .1  denote  Bessel  functions  of  imaginary  argument,  and 

V-Hrp  V"rn*  - ^ J 

h^(t) jh^tt)  arbitrary  functions.  Per  simplicity,  consider  the  special 


ease  v = integer  for  which 


-P 


Then,  it  is  easy  to  pick  out  the  solution  representing  propagation  on 
the  ’outgoing*  characteristics  t+?£  = constant  (note:  with  the  dofinition 
(llj.) , ^ decreases  as  z increases,  hence  the  outgoing  wave  hr s tho  posi- 
tive sign  in  the  characteristic  equation)}  it  is 

(1?)  A =^VL(^p)  j(t), 

where  l(£)  = V(^“) i a^d  j(t)  is  an  arbitrary  function. 

Interpreting  (17),  we  have  the  solution 

(16)  A =Yvj(t+^  - ^)  + - ♦ * 
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where  jn(t)  denotes  the  nth  integral  of  j(t),  the  lower  limits 
being  the  value  of  t + £ at  the  head  of  the  wave.  Hear  the 
head  of  the  wave,  the  successive  terms  decrease  in  magnitude 
provided  ^ 0 , and  the  solution  is  of  exactly  the  same  general 
form  as  (11);  that  is, 

(19)  4><o)f(c)  + <j>(1)f1(£)+...‘, 

where  4 is  the  characteristic  variable  defined  to  be  zero  at  the 
head  of  the  wave , f(4)  is  on  arbitrary  function  (j(t+^)  = f(4)), 
and  the  <j>'s  are  certain  functions  of  the  space  variable.  In  the 
particular  case  (18),  there  are  only  a finite  number  of  terms, 
but  for  non-integral  values  of  v the  series  would  be  infinite. 

We  note,  however,  that  in  this  case  the  scries  is  valid  only 
for  small  4 and  ^ ^ 0*  as**-*  0,  it  is  not  uniformly  valid. 

For  fixed  4,  the  last  term  of  (18)  dominates  for  sufficiently 
small  values  of  ^ $ th-.  order  of  the  terns  should  be  reversed, 
and  the  expression  (18)  (which  in  general  is  replaced  by  an  in- 
finite series)  would  involve  successive  derivatives  of  j (t) . 

Thus,  uc  see  that  although  (19)  is  valid  in  the  initial  stages 
of  the  propagation,  it  may  not  apply  in  later  stages. 

(iii)  Discussion  and  generalization  of  the  results  of  (i)  and  (il) 

An  expansion  >f  the  form  (3.9)  can  be  f ;und  quite  generally 
for  an  equation  of  the  form 

(20)  4»tt  = A2(z)  |<j>Z2  + B(z)  + C ( z ) <|>  j , 

oven  though  the  explicit  solution  is  mt  known,  by  substituting 
(19)  in  the  equation  and  equating  the  coefficients  of  the  f's 
to  zero.  This  was  used  extensively  in  (2j;  in  all  cases  (apart 
from  the  occurrence  of  singularities  as  in  the  problem  discussc-d 
above),  (19)  is  valid  for  small  4>  end  in  certain  cases  such  as 
(11)  its  validity  extends  into  the  region  *v£/r  small.  The  ex- 
pansion is  intimately  connected  with  the  W.K.B.  method  ■£  find- 
ing asymptotic  solutions  of  ordinary  differential  equations. 

For,  the  operational  form  of  (20)  is 

lzz  + B(zHz  +-Jp(z)  - p?/A ^ ( z ) j-  4>  = 0, 
and  the  asymptotic  solution  for  large  p (corresponding  to  small 
4)  is  obtained  in  the  W.K.B.  method  by  the  substitution 


r 


^e‘Pq'z)  ^ -!>(o)(z)  + jk|;(1)(z)  + ...j  f(p), 
where  f(p)  is  an  arbitrary  function  of  p.  It  is  found  that 

r\  

q(z)  = jA“xdz,  hence,  letting  f(p)  be  the  representation  of 
f{t),  the  expansion  is  exactly  (19). 

Another  interesting  feature  of  (19)  is  in  the  significance 
of  Suppose,  there  is  a discontinuity  in  f(^),  then,  since 

f (£)  is  continuous  for  n > I,  (19)  shows  that  an  initial  dis- 
continuity w -'U 1 1 ue  propagated  along  a characteristic  and  its 
magnitude  is  proportional  to  <j/0^  . (Since  <j>  is  the  velocity- 
potential,  discontinuities  of  .-j)  itself  would  not  arise.  How- 
ever,^ is  proportional  to  the  pressure  increase,  and 

|^0^f*(£)  + <j>^^f(£)  + ...  • therefore,  the,  some  result 
is  obtained  f ~>r  pressure  jumps  by  considering  a discontinuity 
in  f'(£))»  0n  the  other  han’,  the  propagation  of  - such  discon- 
tinuities can  be  deduced  by  certain  standard  methods  directly 
from  the  coefficients  of  equation  (20);  if  this  is  J--.no,  exactly 
the  same  function  <j^°^  is  found  f >r  the  magnitude  of  the  dis- 
continuity as  it  moves  along  the  characteristic.  Thus  (19)  al- 
ready contains  the  correct  (linearized)  propagation  of  discon- 
tinuities. Moreover,  since  we  have  a complete  expansion  near  the 
head  .j f the  wave,  additional  information  is  furnished  about  the 
variation  in  the  magnitude  of  th«,  higher  derivatives.  For  ex- 
ample, suppose  as  a special  case  of  (19)  wo  have 

4>  ~ 4(j)(z)  £ + j|'(1)(z)  £2  + ...  . 

Then,  the  liscontinu.iries  in  an-;  <j>  on  £ = 0 are  proportional 
to  and  - A”^c| v ^ j the  discontinuities  in  ^2  nnd  kzz  °n 

E,  = 0 are  proportional  to  d j>  ^ ° Vdz  " n » nnd 

- 2A”1d(})^0  V^z  + A”2.] S'"'  t./dz^  , respe  ctively  and  so  on. 
In  ’epon  ’cntly,  M.  Kline  has  c ; si  bored  a special  case  ;f 
(19)  (an'  -f  the  extensions  which  will  be  -describe  l in  (iv))  * 
he  co:  ..  ).  err;  periodic  waves  an  : thus  f ( £)  = cx  A1  though 
it  is  convi.niijnt  in  our  work  i,  » Iv.k;-  the.  lower  limits  :f  the 
integrals  in  (19)  as  5 *>  0,  the  expansion  is  a solution  ■■£  the 
equation  whatever  the  limits  .are,  or  even  if  the  Integrals  are 
indefinite;  hence,  fn(£)  may  be  taken  as  e1W?/wn.  In  this 
way  Klino  has  an  expansion  valid  for  large  w (for  any  £). 


t 


( i v ) general  problem 
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Thc  extension  of  (19)  tc  include  an  extra  space  variable  is 
now  discussed  and  in  a similar  wav,  any  number  of  space  variables 
could  be  treated.  When  an  additional  space  variable  is  introduced, 
fchoie  no  longer  exists  a union  e sot  of  characteristic  surfaces 
£((t,r,z)  = constant,  since  the  characteristic  conoids  from  any 
curve  ir  the  piano  t = 0 have  an  envelope  which  is  a characteris- 
tic surface.  However,  the  wave  front  t = u(r,z)  is  uniquely  de- 
fined, and  the  obvious  extension  of  the  <g  used  in  (19)  is  £ = 
j?  t - u(r,z),  so  that  at  any  point  in  the  (r,z)  plane,  £ meas- 
ures the  tine  which  lias  elapsed  since  the  passage  of  the  wave 
front.  In  addition,  the  <j>*s  of  (19)  now  become  functions  of  r and 
z.  t\z  before,  this  expansion  will  be  valid  near  the  head  of  the 
wave  since  f /i'n_^  = 0(£)  , and  in  certain  cases  the  region  of  va- 
lidity may  be  extended  to  include  points  at  large  distances,  as 
in  (i).  When  the  polytropic  assumption  is  made  a behaviour  simi- 
lar to  (ii)  is  expected*  in  fact,  near  A=0 , any  wave  will  tend, 
by  refraction,  to  become  parallel  to  the  line  A=0*  hence  the  simi- 
larity to  (ii)  should  be  very  close. 

By  substitution  of  the  extension  of  (19)  into  the  equation 


(8),  wo  have 


(21) 

(22) 

(23) 


ar  + az  = 

2ar  + 2a2  +<°>  + j 

sy.  +‘n>  + a*.  dr> 


i + — + a 
rr  r zz 


££  + a + 
r zz 


= + 
+rr 


f(n-l) 

is + 

r Tzz 


The  first  equation  is  the  well-known  eiconal  equation  which 
must  be  satisfied  by  the  wave  front.  If,  for  example , the  initial 
wave  front  at  t = 0 is  known  i . . the  solution  of  a(r,z)  = 0 is 
given,  (21)  is  sufficient  to  determine  the  function  u(r,z).  Geo- 
metrically, tho  surface  t = a(r,z)  in  the  (t,r,z)  space  is  the 
envelope  of  the  characteristic  conoids  drawn  from  points  on  the 
initial  wave  front,  the  successive  wave  fronts  in  the  (r,z)  plane 
for  different  times  are  the  projections  (on  the  (r,z)  plane)  of 
the  curves  of  intersection  of  the  surface  t = a(r,z)  with  planes 
t constant. 
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It  is  convenient,  in  addition  to  the  set  of  wave  fronts  a = 
constant,  to  consider  the  'rays*  which  are  the  orthogonal  tra- 
jectories of  the  wave  fronts  (pigure  1).  Bre  direction  cosines 
of  a ray  at  any  oint  are  (Aa^,  Aaz)  (by  (21)  the  sum  of  the  squares 
is  1),  and  t.her  .rore  the  rate  of  change  of  any  quantity  along  a 


d . 9 


+ Aa 

z 2 z 


O O 

In  particular,  da/ds  = A(a^+a£)  = l/A,  therefore  the  wave  front  ex- 
pands out  along  tho  rays  with  velocity  A.  In  terms  of  the  ray  deri- 
vative the  left  hand  sides  of  (22)  and  (23)  become  d<}>^n^/ds  + K<j>^n^, 
where  K is  a known  function  of  r and  z,  assuming  that  a(r,z)  has 
already  been  determined.  The  equations  for  tho  t^'s  thus  become 
simply  first  order  rdinary  differential  equations  along  the  rays . 

A.gain,  the  propagation  of  discontinuities  in  the  solution  of 
equation  (8)  can  be  studied  directly  by  an  alternative  method  and 
exactly  equations  (21)  and  (22)  arc  obtained  for  tho  wave  fronts 
and  the  magnitude  of  the  discontiuity . For  this  reason  the  treat- 
ment of  these  equations  - the  theory  of  geometrical  optics  and 
acoustics  (ace  [ Ip 3 ) - is  already  known*  tho  relevant  part  is  repeated 
here.  Equation  (22/  can  be  written 


_i  n_ 


betwocn  any  two  wave  fronts  an1,  two  adjacent  rays,  the  divergence 
theorem  gives 

(25)  /Uar+maz)|(  '}  rRdcT  = C, 

where  the  integration  is  over  the  boundary  of  A BCD,  an'  JL  ,m  are 
the  liroct  :>n  cosines  of  the  ,utw~r!  normal.  Since  (Aa  , Aa  ) are 
the  direction  casinos  of  the  rays,  a +ma  = 0 on  AC  and.  BD, 

x Z 

whole  on  AB  and  CD,  --^a  +ma  = A ( + af  = 1/A.  There f ir*u , (25) 

r z r z 

gives  approx ima t e ly 

2 2 

U(o)  rRAd/A]1-[  l(  /)  rR £<*/h]Q  = 0, 

where  Act  den  >tes  tne  length  of  the  small  segment  of  wave  front 
cut  off  by  the  rays,  and  subscripts  o and  1 indicate  that  the 
quantity  is  evaluated  at  i.B  an::  CD  respectively.  Defining  the 

expansi.n  ro.tio  E = lim  wo  have  rRE/A  is  con- 

Ar  —>o  * 

o 

stant  on  each  ray,  and  in  particular  equals  its  value  at  the  ini- 
tial wave  front.  In  this  way,  wo  may  consider  a(r,z)  and 
<j>^°^(r,z)  to  be  determine:!. 

Equations  (23)  for  the  later  functions  even  though  they 

are  reduced  by  (24.)  to  equati  ,ns  along  the  rays  are  complicated  to 
deal  with.  An  investigation  ■ ! them  v>uld  be  of  interest  in  yield- 
ing additional  information  about  the  region  of  validity  of  (19); 
also,  these  later  terms  show,  in  a similar  way  to  that  described  in 
Oil),  how  dis  continuities  in  higher  derivatives  vary,  according  to 
the  linear  theory. 

II  Improvements  of  the  linearize-1  theory  .and  determination  of 
the  head  shock. 

According  to  the  linearized  theory,  in  the  early  stages  of 
the  propagation  and  even  later  in  certain  cases,  we  have  a solu- 
tion )f  the  form  (19),  and  as  - first  approxinati  ,.n  Wo  take  only 
the  first  term 

(26)  t>  = X(r,z)f  ( t-u(r,z) ) , 

where  f r convenience  in  writing  has  been  replaced  by  X.  It 

is  n w shown  how  (26)  may  be  i.mvr >vcd  and  the  head  sh-clc  deter- 
mined. The  method  will  apply  wh,.mver  -j>  is  of  the  appr  -xim-ate  form 
shown  in  (26)  even  if  it  no  longer  comes  from  ar  expansion  (19).  For 
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examplo,  in  the  problem  d.„jcr  ibed  in  (ii),  (19)  did  not  apply  near 
fy  - o,  bat,  near  = 0,  another  expansion*  whose  dominant  term 
was  3tiil  of  the  f • rm  (26)  (wish  different  X and  f ) , was  valid, 
tbrtefvie,  t \o  method  described  below  could  still  be  applied  in 
that  region. 

The  functions  u(r,z)  and.  X(r,z)  are  assumed  to  be  known  from 
the  geometrical  acoustics  described  in  the  previous  section  and  f 
is  an  arbitrary  function  which  is  fixed  by  applying  some  appropri- 
ate boundary  or  initial  condition.  Unfortunately,  appropriate 
boundary  conditions  would  be  applied  in  regions  outside  the  range 

of  validity  of  the  expansion  (19).  For  example,  if  the  acoustic 

waves  in  1(1)  are  produced  by  the  small  motion  of  a “cylindrical 
piston"  of  initial  radius  equal  to  zero,  the  boundary  condition  is 

applied  where  r is  small,  and  tf  ^^Aor-1(  h(t')dt'  there).  However, 
in  explosion  problems,  the  general  formof  f is  known  from  examina- 
tion of  special  cases  such  as  the  problem  in  [l]j  therefore,  a 
suitable  form  can  be  assumed  when  required. 

Henceforth,  we  shall  be  m.rc  interested  in  the  values  of  u,w, 
and  a than  in  the  velocity  potential  <{>•  these  are  now  set  down. 

'The  dominant  terms  in  u and  w are 

(27)  u = arXF(t~a)  , 

(28)  w = azXF(t-a)  , 

where  F(£)  = -f'(£).  The  sound  speed  a is  obtained  from  Bernoul- 
li’s equation  (4) > but  it  involves  the  pressure-density  relation; 
for  simplicity,  we  assume  that  the  air  is  polytropic,  although  the 
general  ease  could  be  carried  out  in  the  same  way.  Then 

(29)  a A - (J>t  = A + ^ XF(t-a)  . 

In  acc  rdance  with  the  hyp.  ) thesis  which  has  been  developed  and 
discussed  in  detail  in  previous  papers  (Cl],  [2],  [3]),  the  linear- 
ized theory,  which  becomes  inadequate  after  the  disturbance  has 
propagated  for  some  time,  may  be  corrected  quite  simply  to  provides 


For  integral  values  of  v,  the  two  ‘expansions’  were  the  some 
expression  vrrltten  in  different  order*!  from  the  present  point  of 
view,  however,  they  may  be  treated  as  distinct,  since  they  would  be 
in  general . 
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valicl  description  of  the  motion.  In  the  exact  non-linear  theory, 
small  'wavelets’  (whose  paths  arc  characteristic  surfaces  in  the 
(t,r,z)  space)  move  at  the  local  sound  speed  a relative  to  the  fluid 
velocity  (u,w).  The  linear  theory,  however,  approximates  this  value 
by  the  undisturbed  sound  speed.  A,  on  the  grounds  thfet  a-A,  u,  w 
arc  small  compared  to  a.  But , as  a wavelet  travels  the  error  intro- 
duced by  this  assumption  becomes  large,  since  the  small  deviations 
in  the  propagation  velocity  accumulate.  To  correct  this,  the  linear- 
ized characteristic  variable  t-a(r,z)  must  ba  replaced  in  (26)-(29) 
by  (t,r,z),  where  tr  is  a sufficiently  accurate  approximation  to  the 
exact  char act oris t i c variable.  Thus  (2?),  (23),  (29)  become 

(30) 

(31) 


u = u rXF(t), 


w = c^XFfc)  , 

(32)  a = A + Zjr  XP(tt)  , 

where  7g(t,r,z)  remains  to  be  determined . 

A surf  accrzX  t , r , z ) = constant  in  the  (t,r,z)  space  corresponds 
to  th-  propagation  of  a wavelet,  and  gives  a set  of  curves  in  the 
(r,z)  plane,  each  curve  corresponding  to  a single  value  of  t and  rep- 
resenting the  position  of  the  wavelet  at  that  time.  It  is  easy  to 
show  that  the  velocity  of  the  wavelet  normal  to  itself  in  the  (r,z) 

plane  is  a - f,  / ~ The  vol>c.i.  ty  of  the  particles  of  fluid 

t v z ______ 

normal  to  the  wavelet  is  (ut^  + wt^)/' J therefore  the  char- 
acteristic condition,  that  the  velocity  of  the  wavelet  differs  from 
this  particle  velocity  by  the  sound  speed,  gives 


uc  + wtr 
r z 


^Jc2  + ~C2  1/x:2 

v r z Y r z 


± a» 


or 

(33) 


+ utr  + wtr  )2  = a.2(C2  + ts2)  . 


Since  expressions  for  u,w  and  a in  terms  of  tc  are  known,  (33)  "ives 
an  equation  to  determine  tr,  Nov;,  although  TTd iff -..rs  from  t-ct  by  an 
amount  which  becomes  large  (this  being  one  of  the  crucial  points  in 
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tho  discussion),  the  'if foronce  is  still  snail  compared  to  a. 
Hence , we  may  set 


"C  = t - ■-(i’,z)  + p-(r,z)  , 

where  p/~  is  small.  Tin  n,  substituting  in  (33)  and  neglecting  terms 
of  second  order  in  small  quantities,  we  obtain 

(34-) 


an  + 

nr 


ua  + wa  „ . 
u = 2 + Sri 

2^z  . 2 A 


(using  (21)).  From  (30)  , (31)  and  (32),  the  right  hand  side  of  (34.) 
is  ^( Y+l)XFfC)/A^.  Hence,  p.  may  be  taken  equal  to  p(r,z)FCC},  where 

(35)  -rPr  > a28z  ^ 5 


from  (24),  the  left  hand  side  may  bo  replaced  in  terms  of  a deriva- 
tive along  a ray,  to  give 


(36) 


dg.  = X+l  JL. 

ds  2 3 

■i  ± 


To  xfiuo  p uniquely,  it  is  stipulated  that  6=0  on  the  initial 
wave  front)  then  p is  determined  by  integration  along  the  rays.  The 
solutioh  ( 30' » ( 31) » ( 32)  is  n -w  completed  by  the  addition  of  the  im- 
plicit relation  for  "C(  t , r , z ) : 

(37)  t * a(r,z)-|(r+l)P(r,z)FCd)  + TT. 

In  the  problems  of  principal  interest  (in  which  the  disturbance 
is  caused  by  an  explosion)  the  pressure  and  the  sound  speod  increase 
at  the  head  of  the  disturbance)  thus,  F(U)  is  initially  positive. 

It  follows  from  (37)  that  the  characteristic  surfaces  on  which  FfC) 
is  positive  will  eventually  intersect  the  characteristic  surfaces, 
t - a = constant,  of  the  undisturbed  region  ahead  vihcre  FCO  = 0. 
(Another  way  of  expressing  the  same  thing  is  that,  in  the  (r ,7) 
plane,  wavelets  carrying  greater  values  of  Ffc)  travel  faster  and 
overtake  the  others) . If  this  is  allowed  to  occur,  the  solution  is 
no  longer  single  valued  an’  ceases  to  have  physical  significance; 
the  head  shock  must  be  fitto  • in  this  region  in  order  that  the  char- 
acteristic surfaces  meet  it  and  are  cut  off  by  it  bef  re  they  inter- 
sect each  other. 

At  any  time,  the  following  conditions  must  be  satisfied  across 
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the  shock  in  the  ( r,z ) plane: 

(a)  The  component  v.  _ of  th.  particle  velocity  tangential  to  the 
shock  must  be  continuous . 

(b)  v — 2(o.-A)/(y-1)  where  v is  the  component  of  the  particle  ve- 
locity normal  to  the  shock,  v and  a being  measure!  just  behind  tho 
shock. 

(c)  U = i(A+a+v  ) where  U is  the  normal  velocity  of  the  shock. 

' n 2 n'  n 

Those  are  deduced  from  tho  Rankino -Hugoni s t shock  relations, 
neglecting  terms  of  second  and  higher  order  in  the  strength  of  the 
shock.  (?)  and  (c)  are  the  obvious  extensions  of  the  two  shock  con- 
ditions which  were  used  in  [l]  and  (a)  follows  directly  from  the 
conservation  of  the  tangential  component  of  m incntum.  It  may  be 
n-to'  that  (c)  expresses  the  fact  that  the  shock  velocity  is  the 
mean  :f  the  velocity  of  tho  wavelets  on  either  side  of  it;  in  the 
(t,r,z)  space,  it  shows  that  tho  shock  surface  bisects  the  angle  be- 
tween the  characteristic  surfaces  meeting  it  from  either  side. 

It  is  clear  that  the  sh  ck  may  be  assumed  to  be 

(38)  t = a(r,z)  - cr'  (r,z)  sTlr.z). 


where  0/ a is  small,  but  a rain  a crucial  point  is  that  cf  becomes 
largo  compared  to  t-a._  The  component  v^  of  one  particle  velocity 
is  Aua  -Awa  to  the  first  order,  and  this  must  be  zero  by  (a).  But 

(30)  and  (3l)  show  that  the  condition  is  already  satisfied  by  our 
solution.  The  component  v is  Aa  u + Ac.  w to  first  order,  and  from 

(3°),  (3D  this  is  A ( a^+aj" ) XP  fcT)  = XA"1?  ; hence  (32)  shows  that 

I*  z 

(!)  is  also  satisfied.  Tho  final  condition  (c)  determines  the 
shock.  From  (38), 


U = 
n 


+A^(a  cr'  +c.  rf  ) , 
' r r v.  •7.  * 


z z 


since  f/a  is  small;  theref  u'e,  from  (c),  (32)  and  the  above  expres- 
sion for  v , 

(39)  771  • 


In  ad.  lition,  equating  the  values  of  t-a  given  in  (37)  an  : (38),  we 


have 
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d = i3Pfe^  - TT  . 


In  principle,  by  eliminating  (39)  and  (4°)  determine  o' as  a 
function  of  r and  z.  It  is  more  convenient,  however,  to  find  the 
relation  between  r,  z,  amd  the  value  of— cr at  the  shock;  this  re- 
lation specifiesTT as  a function  of  r,  z and  then  (40)  gives 
d\v ,z)  . For  this  purpose,  the  right  hand  side  of  (39)  is  replaced 
from  (35)  by  ^(a„(3_+a  B )FfC)  , and  (40)  is  substituted  for  d to 

L.  i‘  I'  Z Z 

give 

(41)  |(apPr+azPz)F(7r)  + (pF'tc)-l)(ar7^+aTrr)  = 0, 

On  multiplying  by  F(tt)  , (4D  can  be  written 

( r ) 

(U2)  * az  =0. 

Since  the  operator  is  the  derivative  along  the  rays,  the 
term  in  brackets  is  constant  along  each  rayj  but  this  quantity  is 
zero  on  tho  initial  wave  front,  therefore 


p(r,z)  J FfcMdC*  = 0, 


at  all  points  on  the  shock.  It  is  of  interest  to  ooserve  that 
(43)  is  exactly  the  formula  which  was  found  for  the  shock  in  the 
problems  involving  omy  one  space  variable,  except  in  those  cases, 
of  course,  (3  was  a function  only  of'  the  one  variable. 

As  already  mentioned,  (43)  givesTias  a function  of  r and  z 
at  the  shock*  from  it,  otr,z)  may  be  obtained  by  (40 ) and  the 
values  of  u,w  and  a just  behind  the  shock  by  (27),  (28)  and  (29). 
The  strength  of  the  shock  is  found  similarly  from 
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